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Abstract 

The dynamical stability of cosmic rings, or vortons, is investigated for the par- 
ticular equation of state given by the Witten bosonic model. It is found that 
there exists a finite range of the state parameter for which the vorton states 
are actually stable against dynamical perturbations. Inclusion of the electro- 
magnetic self action into the equation of state slightly shrinks the stability 
region but otherwise yields no qualitative difference. If the Witten bosonic 
model represents a good approximation for more realistic string models, then 
the cosmological vorton excess problem can only be solved by assuming either 
that strings are formed at low energy scales or that some quantum instability 
may develop at a sufficient rate. 

PACS numbers: 98.80.Cq, 11.27+d 
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I. INTRODUCTION 



Cosmic strings |1J] have been proposed as seeds for large scale structure formation || and 
as a means to reproduce temperature fluctuations in the cosmic microwave background |3j 
if they appear at the Grand Unified (GUT) phase transition. However, the underlying 
string theories in these calculations assume structureless strings, i.e. the kind for which the 
equation of state is that of Goto-Nambu (having their energy per unit length U equal to their 
tension T, both being strictly constant along the string worldsheet), so the question arises to 
what extent these conclusions are still valid in more complicated cases such as those proposed 
by Witten where a charged condensate is formed in the string's core, thereby inducing 
a current and thus modifying the equation of state into a nondegenerate one. It has been 
argued |S],|5| that in the latter case, equilibrium configurations (called vortons |5[] or rings [Q) 
might exist, that would not radiate all their energy in the form of gravitational waves (which 
is the case for Goto-Nambu strings), and therefore contribute non negligibly to the overall 
matter density in the Universe. Estimates of the corresponding vorton distribution led Davis 
and Shellard || and Carter || to the conclusion that the stability of such states would imply 
a huge matter-density remnant, so that current-carrying strings appearing at the GUT phase 
transition would induce ~ 10 20 times the critical density, and that a premature collapse of 
the Universe would be avoided only if they are lightweight (being produced at energy scales 
at most comparable to that of the electroweak phase transition). 

There are in fact two necessary ingredients for this overdensity problem. The first con- 
cerns the currents in the strings, since this is an essential requirement for the equilibrium 
configurations to exist: rings are centrifugally supported string loops and thus are only de- 
fined when the equation of state is nondegenerate so that there exists a prefered frame in 
which rotation makes sense. This issue has been addressed in a particle physics framework 
with the following conclusion. If the string forming GUT theory is to have as a low energy 
limit the standard electroweak model (a "naturalness" requirement), then the string forming 
fields cannot be arbitrarily decoupled from electromagnetism, and currents will appear in 
strings, either through quantum tunnelling 0], or by spontaneous current generation ||. 
This can in fact be traced back to the nonabelian nature of the GUT model being used: 
since the string forming Higgs field expectation value vanishes in the string's core, some 
charged gauge vectors are massless there, and induce an electromagnetic instability || or 
metastability [|7|, depending on the coupling constant values. In both cases, currents appear 
in strings as a generic feature, so the problem remains. Moreover, it has also been shown 
that even in these realistic (and therefore more complicated) models, the equation of state 
can be well approximated by that of the Witten bosonic model, which was recently calcu- 
lated in detail P|-|TT|]. This is the reason why this particular equation of state is used in the 
present work. 

The second necessary requirement to have a string dominated Universe is that the equi- 
librium configurations be classically stable, at least dynamically. This is not an easily verified 
assumption, and the purpose of this work is precisely to investigate it in the framework of 
the Witten bosonic model, both in the neutral limit [(10[ case where the electromagnetic 
coupling is made to vanish, and in the charge coupled case, by implementing the actual 
value of this coupling constant on the equation of state , including back reaction on the 
constituant fields. Since this work is primarily concerned with dynamical stability, we have 
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not included the electromagnetic corrections either to the the equilibrium condition or to 
the stability constraint, these (expected stabilising) corrections being of a higher order and 
left for further examination [ O . 

This work is arranged as follow. In a first part, we recall the basic dynamical equations 



for a string loop, and the equilibrium condition, as well as the stability constraints ||13 
This is done assuming an underlying equation of state, or, more precisely, by assuming the 
characteristic perturbation velocities to be given. These velocities correspond to transverse 
and soundlike perturbations, being expressible respectively as c 2 -. = T/U and c 2 L = —dT/dU. 
Then, a brief summary of the characteristic features of the Witten bosonic model is given 
and we ultimately apply the stability calculations for this model to show that a finite range 
of currents actually yields stable states. Since it turns out that this range is for low values 
of the current, or for nearly lightlike currents, we argue that most vortons are in fact stable, 
so that, unless quantum radiation provide some efficient way to destabilise them, they can 
only have been created at a low energy phase transition. 



II. STABILITY OF A RING CONFIGURATION 

In order to calculate the dynamical evolution of a cosmic string, we assume it to be 
infinitely thin and characterize the corresponding two-dimensional worldsheet by a set of 
two orthonormal tangent vectors w M and t> M , respectively timelike and spacelike. These 
vectors allow an easy definition of the first and second fundamental tensors of the string's 
worldsheet, namely p^JT5[ 



jf = + K^" = t/V^VatA, (1) 

the latter satisfying the Weingarten identity 

K [H P = 0, (2) 

which is the integrability condition for the existence of a worldsheet containing and 
as tangent vectors. 

Since the string is considered as infinitely thin, its stress-energy tensor T^ v is a distri- 
bution defined on the string worldsheet only, and therefore involves in principle singular 
S— functions. However, use of such distributions can be avoided |15[] by working with a 
surface stress-energy tensor T^, depending on the internal coordinates of the worldsheet, 
r and a say, such that the corresponding distribution valued tensor T^ v can, if needed, be 
obtained from the formula 

T^(x) = {-g)- 1 ' 2 [ dS 2 f^(r, a)5\x - x(r, a)), (3) 



where x(r, a) is a generic point of the worldsheet, dS 2 the surface measure and g is the 
determinant of the metric g^ v of the surrounding four dimensional space-time (which is 
supposed to be flat in the following analysis). Conservation of the stress-energy tensor T^ v 
then implies that the surface stress tensor satisfies a corresponding worldsheet conservation 



law having the form ||I4] 



V\V X T^ = 0. (4) 
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This is the dynamical part of the string motion. It can be further simplified by choosing the 
frame v p ) used in Eq. ([!]) to be identified with the frame in which T is diagonal, with 
eigenvalues U and T respectively the energy per unit length and tension, in the form 

fv> = Uu n u v _ 7W. (5) 

It is convenient for computational purposes to start with and v p as basic independent 
variables whose integration determines the worldsheet as a secondary construct. In such 
an approch |13|], the Weingarten condition (Q) is not satisfied automatically but must be 



included as an additional dynamical equation together with Eq. (f|). Subject to provision of 
the equation of state (see next section), Eqs. (0) and (ffl) provide a complete description of 
the string motion, expressible in terms of the transverse perturbations and the longitudinal 
group velocities 





c T = \ 77 , c L = \ --jrz, (6) 



as the following system: 



K [H P = 0^ V {u<>V p {v v ) - v<>V p {u v )) = 0, (7) 

± liV V p f pv = V KV p (tO - c 2 v"V p {v v )) = 0, (8) 

uJVpTP" = AuPV p (c 2 T ) - (1 - 4K^V P K) = 0, (9) 

v u V p f pu = - AclvPV p (4) + (1 - c 2 T )u uU py p {v v ) = 0, (10) 

where -L fiu = g pu — r\ pv is the orthogonal projection operator and A = — (cf, + c^) -1 . The 



unknown quantities QT6J associated with this system can be chosen as the five independent 
components of the tangent vectors u p and v M and c^, while c 2 L is given a priori as a function 
of <3p by the equation of state. 

Let us now more specifically turn to the case of a circular rotating string with radius R, 
angular speed Q and running velocity v = Rfl, for which one can express the eigenvectors 
in cylindrical coordinates (t, p, 9, z) as 

u" = 7(l,0,«,0) and ^ = 7^,0,1,0), (11) 

where 7 = l/vl — v 2 is the Lorentz factor associated to the running velocity v. It may be 
shown [§,0, that the ring configuration ([□]) automatically satisfies the system of Eqs. (|^)- 



([lOf)) except for Eq. (§) which provides the simple relation (again, not taking into account 
the small long-range electromagnetic back reaction p 



v = c T . (12) 

Configurations for which Eq. (|T2| ) hold are the equilibrium states (vortons or rings) whose 
condition of stability is summarized below. This relation is in fact more general in the sense 



that it must hold also for any non straight equilibrium state with a static Killing vector [JT3 

We now consider the perturbations of the equilibrium state (|T2|) . Without loss of gen- 
erality, and thanks to the symmetries of this configuration, the perturbed quantities can be 
taken in the form of plane waves with pulsation uj and integer angular momentum n, so 
that, in particular, the perturbed parameters describing the string's location read as 
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p = R + SRe i{uJt ~ n0) , (13) 
z = 5ze l{wt - ne) . (14) 

with 5R and 5z constant amplitudes small compared to the unperturbed value R. The 6 
perturbed independent quantities |16] in this particular case are chosen as Su p , Su , Su z , 6u p , 



6u z and out of which one can reconstruct the full perturbed ring, and in particular 5R 
and 5z, this reconstruction being consistent only in the case where the geometrical Eq. (0) 
is satisfied. 

The equations (|7]) - (|T0|) give a linear homogeneous system of six equations with six 
unknowns in which the azimuthal perturbations as exemplified by Su z and 6u z decouple from 
the equatorial ones consisting of the other unknowns 6u p , 6u e , 6u p and fie?. The azimuthal 



part of this system yields the modes |L6 



f 2m; I 

■ £ |^} < 15 > 

which are all stable (their imaginary part vanishing). Therefore, we shall only be concerned 



by the equatorial modes, given [13[] as the solutions of the third degree polynomial in o = to/Q 



v\l + v 2 ){\ - c 2 v 2 )a 3 + 2nv 2 [c 2 - v 2 - 2(1 - c 2 v 2 )]a 2 + 
[Av 2 {\ - c 2 )(n 2 - 1) - (1 + v 2 )(c 2 - v 2 )(n 2 + l)]a + 2n(c 2 - v 2 )(n 2 - 1) = 0, (16) 

with c = cli and where the static mode o = has been implicitely extracted out. This 
equation might have complex roots, in which case the corresponding state will be unstable, 
with characteristic life-time r _1 = |Zmu;|. A previous analysis Hl6| , using explicit solutions 



of Eq. (|IIJ) given by the Cardan formulae, exhibited the stability and instability regions in 
the square (c|, c^) as well as provided analytic expressions for the imaginary part of the 
modes in the latter case. It was found that while approaching the corner c| = c|> = 1 
(where the Witten model is located), the zones where instabilities might develop turn into 
vanishingly thin surfaces so that almost any equation of state must here cross stable zones. 
This observation was the starting point for the following closer examination of the Witten 
model in this context. 



III. THE WITTEN BOSONIC MODEL 

The most simple model that leads to superconducting cosmic strings consists in the 
abelian Higgs model in which a Higgs field $ breaks a U(l) symmetry by means of a non- 
vanishing vacuum expectation value (1$!) = r], thereby giving a mass to a gauge vector 
boson B^, coupled to a simplified representation of electromagnetism in which a charged 
scalar field E (representing some sort of average of all the various possible fields involved in 
the underlying theory) is coupled to the photon A^. The general Lagrangian density that 
describes such fields is 

£ = --\DM 2 - -\D^\ 2 - —H 2 - —F 2 - V($, £), (17) 
2 2 16tt M IQit y ' h K ' 

where the covariant derivatives are defined through 
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= (<9 M + igB^S, ee (0„ + zeA M )S, (18) 

the Maxwell tensors are 

-f^u/ = d^B^, = dfcAv], (19) 

and the interaction potential can be taken as 

v($, e) = ^d$i 2 - v r + m? - v 2 m 2 + + ^\n 2 . m 

This model allows vortex solutions with a S-condensate responsible for a change in the 
equation of state as we now recall. 

In order to study the stability of ring configurations, it is necessary either to calculate 
the full field equations in a circular vortex configuration, or, assuming the thin string ap- 
proximation to hold, to assume the string locally straight for the fields, calculate the various 
integrated (on a transverse plane) components of the stress energy tensor, and plug these 
back into Eq. (|3]). These values can afterwards be used in the previouly developed formalism, 
so it is this approach that we shall follow here. Thus, we shall consider a portion of straight 
string, which we choose to be aligned with a coordinate axis z, and, working in cylindrical 
coordinates, we study a Nielsen-Olesen \TE] vortex solution of the form 



$ = <p(r)exp(iN0), £ = a(r) exp[iip(z,t)], (21) 

for some integer winding number N, ip being possibly restricted to the simple form 

ip(z,t) = az-bt. (22) 

It can be shown [pTCj , pTlJ] that the only arbitrary parameter for a string in this model (assuming 
the underlying coupling constants to be given) is the phase gradient 



v 



±\ZlWl, (23) 



where the + or — sign must be chosen according to whether the conserved electromagnetic 
current 

J" ee = a 2 (V^ + eA») (24) 

e dAf, 

is respectively spacelike or timelike. A very important point concerning this current, and the 
corresponding existence of centrifugally supported equilibrium states, is that it is defined 
even in the so-called neutral limit (TtJ where the electromagnetic coupling constant e is 
made to vanish. This model therefore accounts in particular for neutral-carrier condensates 
in cosmic strings, and allows an easy recognition of the fact that the long range feature of 
electromagnetism is essentially irrelevant in ring dynamics. On the figures explained below, 
curves have been plotted using both e = and a large value for this parameter producing 
only minor quantitative corrections (exagerated on the figures). 

The other conserved quantity that is needed is the stress energy tensor 
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(25) 



which yields, in the notation of the previous section, Eq. (|) 

U = f tt = 2ir J rdrT a (26) 

and 

T = -f zz = -2tt f rdrT zz . (27) 



The procedure now goes as follows: to each value of the state parameter v corresponds a 
(numerically computed) field configuration which can be integrated to yield U and T, with 
which we calculate the velocities ct and c^. (More details concerning these computations 



can be found in particular in Refs. |I0,P|.) A characteristic result is shown on Fig. 1 on 
which are plotted these velocities as functions of the state parameter v expressed in units 
of the mass of the current-carrier S, in the neutral limit (dashed curves) and in the case 
where the full electromagnetic self action has been taken into account (full curves). Again, 
it should be stressed that the corrections resulting from the inclusion of the electromagnetic 
coupling constant e affect our results only in a quantitave way, the perturbation velocity 
plots being almost indistinguishable in most of the parameter space. 



IV. RESULTS AND CONCLUSIONS 

The Witten bosonic model has several underlying parameters, all of which are supposed 
to be fixed by the underlying string-forming microscopic theory. For each particular string, 
there is only one that remains independent, namely the squared phase gradient v 2 of the 
current- carrier condensate, called the state parameter, whose sign reflects the timelike or 
spacelike nature of the superconducting current, and whose amplitude gives, in a nontrivial 
way, the amplitude of the corresponding current, and the degeneracy of the stress-energy 
tensor. For various values of the underlying parameters, we have derived the variations 
of the energy per unit length U and the tension T with the state parameter, enabling us 
to calculate the actual values, in these models, of the perturbation velocities ct 2 = T/U 
and cl 2 = —dT/dU, as illustrated on Fig. 1, both for vanishing electromagnetic coupling 
constant and with inclusion of the full self action on the string-forming fields. We have 
then used these values to question the stability of ringlike configurations (vortons) against 
azimuthal and equatorial perturbations, and we found that generic results, meaning ones 
roughly independently of the values of the underlying parameters (coupling constants A^, 
X a , q, e, masses, • ■ ■), could be drawn as exemplified in Figs. 1 and 2. 

The first point, as was already emphasized in previous work (but whose implications had 
not yet been worked out in full detail), is that the velocity of transverse perturbation always 
exceeds that of longitudinal perturbations, even for very low current values, so that the first 
order approximation for the evolution of a superconducting string, namely that for which 
an effective action of the form 

S = const f d 2 s(y^ip + eA^) 2 , (28) 
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where the transverse degrees of freedom have been integrated out, is in fact unsatisfactory 
as soon as one needs to consider derivatives of the equation of state. For such purposes, the 



alternative action 15 



S = const J d 2 s^Jm 2 + (V^ip + eA M ) 2 , (29) 

with m a constant with the dimension of mass, should preferentially be used since it re- 
produces most of the actual features of the Witten model, and besides, is completely in- 
tegrable (in the sense that analytic solutions for the string worldsheet can be explicitly 
constructed (19|.|2(|) and thus leads to much simpler equations of motion. This may be un- 
derstood by stating that the crucial contribution for the derivatives involved in c| is given 
by the fourth order term in the phase gradient, assumed to be identically zero in the ac- 
tion (p8~D, though not in Eq. (|2~9|) . More precisely, it can be seen that Eq. ( p8]) implies cl = 1, 
whereas in fact, in the Witten model, it is not only less that the speed of light, but also less 
than ct- In that sense, Eq. fl2"9|) gives the better approximation cl = op, although stability 
considerations do not apply since this relation has been shown to imply absolute stabil- 



ity against any dynamical perturbations. For currents small enough, and when the internal 
degrees of freedom of the strings are neglected, the approximation (p8|) can still be used for 
some purposes, but only in a more restricted range of applications than previously thought. 

Now the problem is that if the longitudinal velocity had been greater than the transverse 
one, the ring dynamical stability issue could have been addressed far more easily since in 
this region of the (cl,ct) plane, circular configurations are always stable. However, it is 
less obvious what will happen in the Witten string case: the regions of instability in the 
(cl,ct) plane reduce to lines close to the cl = or = 1 edge where the most interesting 
part of the Witten-model's equation of state is. Therefore, although it is not possible to 
show it explicitely on a plot (the number of unstable zones is in principle infinite as one 
approaches the cy = 1 line and it was not possible to draw a clear graph), it appears 
that the equation of state will cross various stable and unstable regions as it goes away 
from the stable point cl = ct = 1 (including in particular the null current case), before 
ultimately reaching the unstable region. This can be shown alternatively as on Fig. 2 where 
the characteristic inverse life-time of the configuration, expressed in units of the angular 
velocity Q of the ring, is calculated [as the larger imaginary part of the solutions of Eq. (|T5D] 
as a function of the state parameter. This figure shows explicitely the range of phase gradient 
where the corresponding ring states are stable. Various remarks need to be made. First, the 
electromagnetic correction to the equation of state yields a small correction to the actual size 
of the stability region, but does not modify the result qualitatively. It should be emphasized 
that for a realistic underlying field theory in which the mass of the string forming Higgs 
particle is expected to be much larger than m a , the relevant coupling constant would be so 
small |TJgg 

e 2 (m ff /M ) 2 ~ 1(T 6 (30) 

that its effect would be imperceptible. In order to obtain an effect large enough to be 
visible on the figure, an artificially exagerated value e 2 (m a /M ( f ) ) 2 = 0.1 has therefore been 
used. Therefore, one important conclusion that can again be drawn from these results is 
that electromagnetic corrections to the field equations in the string core do not significantly 
modify the actual dynamics of a string, so that use of the neutral limit model |l(J is justified. 



S 



The last point we want to emphasize concerns the vorton excess problem. As the equi- 
librium condition happens to be reached primarily for quasinull currents [||], for which v is 
very small, it can be argued, on the basis of our results, that most of these vortons are in fact 
dynamically stable since this is precisely the region where stability actually occurs. So our 
final conclusion is that if the Witten bosonic model is a good approximation for describing 
superconducting cosmic strings, then any string-forming theory will produce dynamically 
stable vortons which can lead to an overdensity in the Universe unless the phase transition 
happens at low enough energy (estimated of order 10 TeV 0), or if quantum instabilities 
can develop with sufficient rate. This can also be understood as a constraint on the scheme 
of symmetry breaking in any string forming GUT model or on the values of its microscopic 
parameters. 
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FIGURES 



FIG. 1. Plot of the transverse and longitudinal velocities = T/U and c| = —dT/dU as 
functions of the phase gradient v 2 = Id^tpd^tpl of the current carrier (in units of the inverse cur- 
rent-carrier mass). Dashed lines are for the neutral limit with vanishing electromagnetic coupling 
constant, whereas solid lines include the full electromagnetic back reaction on the underlying fields, 
assuming an unreasonably large value for the coupling constant e to enhance the effect. 

FIG. 2. The inverse life-time r _1 = |Xma>| given by the complex roots of Eq. (16) for the Witten 
bosonic model, in units of the angular velocity Q of the ring. A finite region of the state parameter 
space has an infinite life-time (i.e. yield stable rings), for small or quasi-lightlike currents, and this 
is a generic feature in this model. 
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